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ABSTRACT 

We describe an algorithm that measures self-consistently the relative galaxy contribu¬ 
tion in a sample of radio-quasars from their optical spectra alone. This is based on a 
spectral fitting method which uses the size of the characteristic 4000A feature of el¬ 
liptical galaxy SEDs. We apply this method to the Parkes Half-Jansky Flat Spectrum 
sample of Drinkwater et al. (1997) to determine whether emission from the host galaxy 
can significantly contribute to the very red optical-to-near-infrared colours observed. 
We find that at around 2cr confidence, most of the reddening in unresolved (mostly 
quasar-like) sources is unlikely to be due to contamination by a red stellar component. 
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1 INTRODUCTION 

The recent discovery of large numbers of radio-selected 
quasars with very red optical-to-near-infrared colours sug¬ 
gests that existing quasar surveys may be severely incom¬ 
plete. Based on a high identification rate in the optical and 
near-infrared, Webster et al. (1995) found a broad range 
of colours with 2 < R — A < 10 for fiat-spectrum ra¬ 
dio quasars in a subsample of the Parkes 2.7GHz survey 
(Drinkwater et al. 1997; hereafter ‘Parkes quasars’). For 
comparison, quasars selected by standard optical techniques 
show a small scatter around B — K ^ 2.5. Two theories have 
been proposed to explain the large scatter: first, Webster et 
al. (1995) interpreted this in terms of extinction by line- 
of-sight dust (see also Masci 1997), and second, Serjeant & 
Rawlings (1996) suggested that this was due to intrinsically 
red optical/infrared synchrotron radiation associated with 
the radio emission. This paper explores a third possibility: 
that the red colours are due to contamination by starlight 
from the host galaxies. 

In a recent near-infrared imaging study of a sample 
of radio-quasars selected from the low-frequency (408MHz) 
catalog, Benn et al. (1998) found that sources with red B — K 
colours to z ~ 2 could be explained by an excess of host 
galaxy light in K. Most of the images appeared non-stellar 
(or resolved) suggesting that indeed starlight was responsi¬ 
ble for the redness in B — K colour. It is important to note 


however that all these sources were associated with extended 
radio galaxies of which a majority are of the steep-spectrum 
type, a common feature of low-frequency selected samples. A 
large fraction are also often associated with luminous giant 
ellipticals. In view of the simple orientation-based unified 
model for radio-loud AGN, it is possible that galaxy light 
is more easily detected in these sources due to anisotropic 
obscuration and/or unbeamed emission of the central AGN. 
The radio-quasars explored in this paper are all of the flat- 
spectrum type selected at moderately high frequency. Thus, 
they are not expected to exhibit similar properties in the 
near-infrared. 

The B — K colours of normal radio galaxies are known 
to be quite red, exhibiting a similar dispersion to those ob¬ 
served in the Parkes sample. These sources are often associ¬ 
ated with giant ellipticals and their colours appear to be uni- 
forml y distributed over the range 3<R — A< 7 for redshifts 
2<2 ( McCarthy 1993 and references therein| ). It is possible 
that the host galaxies of Parkes quasars also exhibit simi¬ 
lar properties. To determine the importance of host galaxy 
light in Parkes quasars, we need to quantify its contribution 
throughout the optical to near-IR wavelength region. 

Determining the host galaxy properties of radio-quasars 
and BL-Lac-type sources are crucial for studies of AGN 
evolution and testing unified schemes. A significant num¬ 
ber of BL-Lacs and other core-dominated radio loud AGN 
are surrounded by nebulosities whose o ptical spectra are 
very similar to those of giant ellipticals (^g. Miller, French 
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fc Hawley 1978|; Ulrich 1988; Taylor et al. 1996). A ma¬ 


jority of these host galaxies have been detected in sources 
at relatively low redshifts, 2 ^ 0 . 1 , where the galaxy flux is 
easily detected. Since most quasars are at high redshifts. 
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host galaxy detections have been difficult due to contami¬ 
nation by their strong nuclear emission. Recent high resolu¬ 
tion imaging using HST has revealed that a ma jority of high 
redshif t quasars also reside in gia nt ellipticals (Hutchings & 
Morris 1995; Bahcall et al. 1995). 

The optical spectral energy distributions (SEDs) of el¬ 
lipticals are all similar in shape. At A^SOOOA, they rise 
steeply up to ~ 4000A where a sharp step-like cutoff is ob¬ 
served. This cutoff is often refered to as the “4000 Angstrom 
break” and is caused by the dearth of hot and young, usu¬ 
ally O and B-type stars with time, whose spectra mostly 
dominate at wavelengths from « 4000 A into the UV. For 
A^4000A, the existence of a multiplet of heavy metal ab¬ 
sorption lines arising from stellar atmospheres, results in a 
steepening of the break. The presence of this feature thus 
provides a signature for determining whether an underlying 
elliptical host is contributing to the total fight. 

In this paper, we describe a new unbiased method that 
uses the size of this characteristic break in the optical spectra 
of Parkes quasars to quantify the host galaxy contribution. 
We use this method to investigate whether the host galaxies 
of Parkes quasars can significantly contribute to the optical 
and near-IR continua and hence cause the spread in B — K 
colours observed. 

This paper is organised as follows: In section we out¬ 
line our method used to quantify the host galaxy contribu¬ 
tion. A discussion of the spectral data and assumed input 
parameters is given in section Results are presented in sec¬ 
tion Implications for AGN unified schemes are explored 
in section ^ and results discussed in section All results are 
summarised in section |^. 


2 OUTLINE OF METHOD 

The galaxy contribution can be estimated using an algo¬ 
rithm that determines the relative strength of the character¬ 
istic 4000A break in the optical spectra of Parkes sources in 
an unbiased way. No other feature in a generic quasar spec¬ 
trum will mimic this feature. In this section, we describe this 
algorithm. 


2.1 Assumptions 



A/Angstroms 


Figure 1. Composite optical spectrum of 18 Parkes sources which 
appear spatially extended in Bj (solid curve), and a model spec¬ 
trum for an ‘old’ elliptical SED from Bruzual & Chariot (1993) 
(dashed curve). 



A/Angstroms 


Figure 2. Spectral evolution of the IGyr burst model of 
Bruzual & Chariot (1993). The age in Gyr and the approxi¬ 
mate corresponding redshift for a formation epoch Zf = 5 with 
{qo, /i5o)=(0.5,l) is shown beside each spectrum. 


There are two important assumptions that will be required 
by our algorithm. First, we assume that the spectral shape 
defining an underlying quasar optical continuum is ‘smooth’ 
and contains no breaks. This is justified by our current 
knowledge of quasar optical spectra. Our choice for this 
shape is somewhat arbitrary, and will be further discussed 
in section ^.2| . 

Second, we need to assume the input galaxy spectrum 
which defines the shape of the 4000A break to be used in 
our analysis. We have formed a composite optical spectrum 
of 18 Parkes sources all of which appear ‘spatially extended’ 
in Bj. This composite is shown as the solid curve in Fig. ^ 
and clearly shows the characteristic 4000A break. Superim¬ 
posed (dashed) is an elliptical galaxy SED predicted from 
the stella r population synthesis models of Bruzual & Char¬ 
iot (1993). This model is for a galaxy of age 8Gyr, and as¬ 
sumes a IGvr burst described by a Salpeter initial mass 
functio: i with no star formation thereafter. In fact, any ‘old 


aged’ model with long completed star formation will in gen¬ 
eral be similar. Models assuming ‘long lived’ constant star 
formation rates lead to considerably bluer optical-UV con¬ 
tinua, characteristic of those observed in spirals. The com¬ 
posite and model in Fig. are in excellent agreement, except 
for the emission fines arising from ISM gas in the composite. 
In our analysis, we shall assume that this model represents 
the underlying elliptical SED in each of our Parkes sources. 

For sources with redshifts z>l, the 4000 A break is red- 
shifted out of the wavelength ran ge available in the spectra 
of Parkes quasars (see section 3.1) and thus, cannot be used 


in our algorithm. We are therefore restricted to a < 1. For 
simplicity, we assume that the general shape of our model el¬ 
liptical SED is independent of redshift to 2 ~ 1. The colours 
of radio galaxies observe d to 2 ~ 3 are consistent w ith forma¬ 


tion redshi fts 4<2/<20 (Spinrad & Djorgovski 1987; Dunlop 


et al. 1989 ). These formation redshifts are also suggested by 
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models of gala xy formation in the CDM scenario (3g. White 
& Frene 1991). The spectral synthesis models predict that 
the general form of an elliptical SED shown in Fig. ^ can be 
immediately established following an almost instantaneous 
(~ IGyr) burst of star forma tion. Adopting the IGyr burst 
model of Bruzual & Chariot (1993) (updated 1995 models), 
we show in Fig. ^ the spectral evolution as a function of 
age since the initial starburst. Also shown beside each age 
is the approximate redshift for a formation epoch of Zf = 5. 
Initially, the 4000A break region is diluted by starburst ac¬ 
tivity. After a few gigayears, the stellar population evolves 
passively, maintaining an almost uniform SED shape (see 
Bruzual & Chariot 1993). For formation epochs z>5, it is 
apparent that the generic elliptical SED shape of Fig. ^ is 
easily established by 2 ~ 1. 


2.2 The Algorithm 

In general, the total flux at a given wavelength observed 
in a Parkes source, /t(A), can be modelled as the sum of 
light contributed by the central quasar or AGN, fq{X), and 
any underlying host galaxy fg{\) (eg. Fig. |^. We write the 
relationship between these quantities as: 


MX) = fT{X)-Cfg{X), 


( 1 ) 


where c is a scaling factor giving an arbitrary measure of the 
amount of galaxy light we wish to determine for a particular 
source. 

Given /t(A) for a particular source and an arbitrary 
galaxy spectrum /g(A) (Fig. |^), our aim is to determine the 
value of c such that fq{X) looks something like a quasar 
spectrum. From our knowledge of quasar optical spectra, 
an obvious choice is to require that fq{X) be “smooth” 
and contain no breaks. Thus, the basis of this algorithm 
involves subtracting an arbitrary amount of galaxy flux, 
cfg{X) from /t(A) such that the resulting spectrum fq{X) 
appears smooth (see below). When this is achieved, the frac¬ 
tion of total light at a given wavelength contributed by the 
host galaxy can be estimated by normalising: 

In order to implement the above algorithm, we need 
to define an acceptable form for the shape of the quasar 
spectrum fq{X). Our only requirement is that this spectrum 
be smooth and hence our choice is somewhat arbitrary. We 
choose fq{X) to be a power-law, parameterised by: 

/.(A) = /g(A^..)f-^)“, (3) 

where the slope a. is determined between two fixed wave¬ 
lengths, Xmin and Xmax (see below): 


In \^fq IXmin )//.(A max')\ 


In [Amin/A max] 

Using Eqn. hi, the fluxes defined in Eqn. W can be written: 


( 4 ) 


/.(A min ) = /t(A min )-c/g(A min ) 
/.(A max ) = MX max )-Cfq{X max ) • 


( 5 ) 


For a discussion on how the wavelengths Amin and Amax axe 
chosen and the fluxes in Eqn. ^ measured, see section 3.1 



Wavelength 


Figure 3. Definition of our “figure of merit” A. For a particular 
source spectrum /t(A) (whose 4000A break region is shown ex¬ 
agerated), we subtract an amount of galaxy flux, fg{X) (Fig. ^ 
until A is a minimum and a smooth power-law fq{X) results. 


The validity of our assumption of a single power-law for 
fq{X) is discussed in section ^ 

To apply this algorithm in a self consistent way to each 
of our optical spectra, we need to define a figure of merit in¬ 
dicating the point at which the maximum amount of galaxy 
spectrum has been subtracted and a “smooth” quasar spec¬ 
trum (ie. a power-law) is achieved. Let us first consider the 
rest frame optical spectrum of a source, /t(A), suspected 
of containing a 4000A break. This is illustrated in Fig. ^ 
Furthermore, let us consider a “smooth” power-law (ie. the 
underlying quasar spectrum fq{X)) between two wavelengths 
Amin and Amnx ou either side of the 4000A break region as 
shown. We define our figure of merit as representing the 
area A of the shaded region in Fig. As certain amounts 
of galaxy spectrum, cfg(X) (where fg{X) is given by Fig. 
are gradually subtracted from /t(A), A will decrease and be¬ 
comes a minimum when the break disappears. We can thus 
determine the value of c when this occurs, allowing us to 
estimate the fractional galaxy contribution from Eqn. 

For a given amount of subtracted galaxy flux, c/g(A), 
we can write A in terms of c from Eqn. ^ as follows: 

A = ^ '”“*|/g(A)-/T(A) + c/g(A)| dX. (6) 


With fq (A) defined by Eqns. H, M and m A can be written: 


A = 


[ IfriXmax) - cfg{Xmax)] 

J \ ^max ' 




- MX) + C/g(A) 


dA, 


( 7 ) 


where a{c,fT,fg) is defined by Eqns. ^ and ^ Thus, one 
would only have to minimise A with respect to c in order to 
determine the maximal galactic contribution. We introduce 
however, an additional factor in Eqn. ^ whose purpose will 
be to make best use of the available data and optimise our 
algorithm. 
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2.2.1 Data Optimisation 

Due to systematic effects, each Parkes optical spectrum suf¬ 
fers from a considerable amount of noise at the range of 
observed wavelengths. Although the noise may have some 
wavelength dependence, we assume it is constant. Given 
this assumption, we have used an optimal method that gives 
more weight to those wavelengths in a particular spectrum 
where the residual |/g(A) — /t(A)| (see Fig. is likely to be 
a maximum. In other words, the galaxy subtraction process 
will depend most sensitively on the spectral shape around 
the 4000A break region. Our assumption that the noise is 
wavelength independent will greatly simplify our optimisa¬ 
tion method, since otherwise, we would have to simultane¬ 
ously optimise those observed spectral regions with the high¬ 
est signal-to-noise. 

The basis of this method simply involves convolving the 
integrand in Eqn. ^ with a weighting function G(A) that 
gives more weight to regions on either side of the 4000A 
break within Xmin < A < Xmax (see Fig. ^). We choose to 
define G{X) imrely from the galaxy spectrum /g(A) (dashed 
curve in Fig. W). This is defined as the difference (or residual) 
between a power-law and the galaxy spectrum within the 
range Xmin < X < Xmax, containing the 4000A break: 

G(A) = |/gPL(A)-/g(A)|, (8) 

where 

foPLiX) = /g(A™.)(-^)““ 

\ ^max J 

and 

)//3(A max)] 

In [Amin/A max] 

The function G(A) will peak at wavelengths on either side 
of the 4000A break. Thus, by convolving G(A) with Eqn. ^ 
relatively more weight will be given to spectral data at 
these wavelengths, where our algorithm mostly sensitively 
depends. 

Thus, the function to minimise can now be written: 


A = 


r^max / \ \ 

/ [hiXmax) - cfg{Xmax)] (v- ) 

J \ ' ^max / 


\ -^aOjT.fg) 


- MX) + C/g(A) 


G(A) dX, (9) 


where G(A) is defined by Eqn.|^. Our aim is to determine the 
value of c for a particular source spectrum /t(A) such that 
A is a minimum. In FigAd, we show examples of the figure of 
merit function A (Eqn. ^ as a function of the parameter c for 
two sources whose spectra are shown in Eig. If the weight¬ 
ing factor G(A) in Eqn. ^ is neglected, we find that the values 
of c at which A is a minimum decrease by up to 5% in most 
sources where the break feature appears relatively weak or 
absent. In sources with strong 4000A breaks, there is no 
significant difference. Given the value of c that minimises A, 
the fractional galaxy contribution at some wavelength can 
now be computed using Eqn. 


2.3 Error Determination 


In spectra where the galaxy contribution is relatively weak, 
an upper limit on its contribution at some appropriate level 
of significance would be required. In order to do so, both 
random and systematic errors need to be investigated. The 
magnitude of these two types of errors are estimated and 
compared in section 4.2 Here, we briefly outline the method 
in their determination. 


2.3.1 Random (Statistical) Errors 

The random error at some confidence level in the galaxy 
contribution is estimated by computing the value of c corre¬ 
sponding to the statistical error in A (Eqn.|^. Since Eqn. ^ is 
actually a discrete sum over wavelength bins Ai, from Xmin 
to Xmax, the statistical error is determined by adding the 
error for each individual bin in quadrature, so that 


a(A) = 






( 10 ) 




where (Ji(I) is the error in the integrand / of Eqn. y for bin 
Ai. (Ji{I) will depend only on uncertainties in the measured 
source fluxes friXmax) and fT{Xmin) in Eqn. y. As shall 
be discussed in section 3.1, these fluxes are estimated by 


calculating the median continuum flux in wavelength bins 
centered on Xmax and Xmin- We estimate the corresponding 
uncertainties by computing the rms deviation from the mean 
spectral flux in these wavelength bins. 


2.3.2 Systematic Errors 


Our definition of the underlying quasar continuum fq{X), in 
each source only requires that it be smooth and contain no 
breaks. This implies that the shape of fq{X) is somewhat 
arbitrary and thus it is possible that our quantitative mea¬ 
sures of the galaxy contribution may strongly depend on its 
assumption in our algorithm. The uncertainty introduced by 
this possible systematic effect will be investigated. 

All our calculations assume that fq{X) is a power-law 
(ie. Eqn. 1^), since the continua of optical quasar spectra are 
well represented by fn oc where a ^ 0.2 — 0.3 (Francis 


199C). To investigate the effects of assuming a different form 


for fq{X) however, we also apply our algorithm by assuming 
for simplicity that fq{X) is a straight line joining Xmin and 
Xmax in Fig. M. This is parameterised as follows: 


fq(.X)line — 


[/.(A max ) fq{Xmin)] 

Xmax Xmin 


(X-X max ) + /^(A max)? (11) 


where fq{Xmin) and fq{Xmax) are defined by Eqn. g. A com¬ 
parison in the fractional gala^ contributions resulting from 
our use of a power-law (Eqn. ^ and a straight line (Eqn. [ll| ) 
for fq{X) in our algorithm, will enable us to estimate the 
magnitude of this systematic effect. Results are presented in 
section 4.2. Other possible sources of systematic error are 


also discussed in this section. 


2.4 Summary 

To summarise, we have presented in this section a method 
to determine the relative galaxy contribution in each Parkes 
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Figure 4. Examples of our figure of merit function A (defined by Eqn. g) as a function of the parameter c. Source spectra are shown 
in Fig. |. 


source in a robust way. Our algorithm requires the follow¬ 
ing two input assumptions: First, an elliptical optical SED, 
/g(A), defining the shape of the characteristic 4000A break. 
This we assume is a model from Bruzual & Chariot (1993) 
(see Fig. |^. Second, we require a spectral shape defining the 
underlying quasar spectrum fq{X)- Unless otherwise speci¬ 
fied, all our calculations shall assume a power-law for /q(A). 

With the only requirement that fq{X) be smooth and 
contain no breaks, the “suspected” 4000A break in each 
source spectrum /t(A), is subtracted until the residual be¬ 
tween fq{X) and galaxy subtracted source spectrum is a 
minimum. The galaxy contribution is estimated from the 
amount of galaxy, /g(A) subtracted. To apply this proce¬ 
dure in a self-consistent manner, we have defined a figure of 
merit given by Eqn.|^. This is minimised with respect to the 
parameter c from which the fractional galaxy contribution 
can be easily computed using Eqn. 


3 DATA AND INPUT PARAMETERS 


Out of the 323 sources in the Drinkwater et al. (1997) sam¬ 


ple, we have optical spectra for 194 or 60% of the sam¬ 
ple. Some of these are from the compilation of Wilkes et 
al. (1983) and others are from recent observations on the 


AAT and ANU 2.3m (see Drinkwater et al. 1997). For these 
latter observations (65 sources), the spectra cover the ob¬ 
served wavelength range: 3200A<A<10000A with a resolu¬ 
tion of ^ BA in the blue and ~ 20A in the red (>5200 A). 
Other spectra have typically a mean resolution ~ lOA and 
cover the range: 3200A<A<8000A. 


3.1 Applying the Algorithm 

Before applying our algorithm, each Parkes optical spectrum 
is redshifted to its rest frame. We then define the fixed rest 
frame wavelengths Xmin and Xmax within which the sus¬ 
pected 4000A break and our figure of merit (Eqn. ^ is de¬ 
fined (see Fig. This wavelength region is chosen using 
the following criteria: first, the highest redshift sources will 


have the rest wavelength Xmax redshifted out of the obser¬ 
vational wavelength range of the spectra. These sources will 
not be able to be used in our algorithm. We therefore need 
to choose Xmax such that the number of sources in which 
our algorithm can be applied is not significantly reduced. 
Second, we need a wavelength range Xmin < A < Xmax that 
makes ‘optimal’ use of the shape of the 4000A break region 
defining our figure of merit A (the residual area in Fig. 

In other words, we need to ensure that this region is un¬ 
ambiguously defined and clearly represented in each source 
spectrum. As a compromise, we assume Xmin = 3500A and 
= 5080A in every source. 

Relative measures of the fluxes: fTiXmin,max) and 
fg{Xmin,max) in the source and galaxy rest frame spectra 
respectively (see Eqn. are determined as follows. We first 
define wavelength bins of width ~ 200A and ~ 400A cen¬ 
tered on Xmin and Xmax respectively, and then calculate 
the median continuum flux in each bin. The wavelength re¬ 
gions defining these bins however may contain absorption 
and emission lines. Such lines are likely to bias our esti¬ 
mates of the continuum level in these regions. From the bin 
widths defined above, the rest wavelength regions of inter¬ 
est are: 3400A < A < 3600A and 4880A < A < 5280A . 
From the available source spectra, we find that no lines 
are likely to contaminate the short wavelength bin. For the 
long wavelength bin however, we find that the emission line 
doublet [OIII]AA4959,5007 and a weak absorption feature 
at ^ 5170A (possibly from Mgl) are strong contaminants. 
To avoid significant contamination, our algorithm excludes 
regions of width lOA centered on these lines. 

Given the definitions above, the rest wavelength range 
required by our algorithm will be 3400A<Arest;$5280A. 
With a maximum observed wavelength of Amaa:(obs) ~ 
lOOOOA in ~ 30% of the available spectra, we find that only 
sources with redshifts ^<0.9 can be used in our algorithm. 
In total, we have about 53 spectra in which our algorithm 
can be applied. 


4 RESULTS 
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4.1 Spectral Fits 

With the input parameters from the previous section, Eqn. ^ 
is minimised numerically with respect to the parameter c for 
each spectrum. Having found the value Cmin that minimises 
Eqn.^, we can reconstruct the initial source spectrum /t(A) 
around the 4000A break region within 3500A(Amin) < 
A < 5080A(Amaa:). This is done by superimposing the max¬ 
imum amount of galaxy spectrum Cminfg{\) generated by 
the algorithm, and a power-law representing the underlying 
“smooth” quasar spectrum fq{X). From Eqns. ^ and ^ these 
reconstructed model spectra can be represented: 


/t(A)„ 


— CminfgiX) -j- fq{X) 

= Cmin FW + [MX 


/ X 
V ^max J 


max) Cmin 
JT,fg) 


/s(A maa:)] 


( 12 ) 


A visual comparison between these model and observed 
spectra in the range Xmin < X < Xmax will allow us to in¬ 
vestigate the accuracy of our algorithm in reproducing the 
observed spectra around the 4000A break region. Recon¬ 
structed spectra fT{X)modei are compared with the observed 
spectra, /t(A), for a number of sources in Fig. ^ On each 
spectrum, we also show our power-law fit fq{X). Since the 
power-law spectral index for fq{X) is defined purely within 
Xmin < X < Xmax (as required by our algorithm; see Eqns ^ 
and 1^, a model spectrum will attempt to accurately re¬ 
produce that observed region which falls in this wavelength 
range only. As expected, observed spectra showing strong 
4000A breaks with SEDs similar to that given in Fig. are 
reproduced very accurately about the break region. 

In a majority of observed spectra where no signihcantly 
strong breaks are discernable to the human eye however, our 
algorithm nevertheless attempts to ht for a break. Unbiased 
estimates of the galaxy contribution using the relative s izes 
of these breaks, however weak, are presented in section 4.3 


4.2 Systematic versus Random Errors 


As discussed in section 2.S, a possible source of systematic 
uncertainty is in our assumption of the shape of the un¬ 
derlying smooth quasar continuum fq{X). To explore this, 
we compare the relative galaxy contribution obtained by 
assuming hrst, a power-law (PL) (Eqn. for fq(X) and sec¬ 
ond, a straight line (L) (Eqn. |ll[ ). Using our algorithm and 
these two definitions for fq{X), we have computed the frac¬ 
tional galaxy contribution at 5000A (rest frame). Results 
are shown in Fig. ^ 

At first glance, estimates for the galaxy fraction using 
the power-law and straight line for fq{X) agree very well. 
There is relatively little scatter about the diagonal line defin¬ 
ing the equality fracr, = fracpi, except for a distinct popu¬ 
lation with fraci:, < fracpL. No distinguishing feature in the 
optical spectra of this latter class is immediately apparent. 
It is likely that a power-law (rather than a straight line) 
within Xmin < X < Xmax for these sources provides a better 
representation of our figure of merit A in Fig. ^ We quantify 
the systematic error from the RMS scatter in the difference: 
S — fraci — fracpi, which we denote by (^((I). For this sys¬ 
tematic effect, we therefore estimate a la uncertainty in the 
galaxy fraction at 5000A of at most afmc. — 0.02. 



gal frac.(power law) 


Figure 6. Fractional galaxy contribution at SOOOA (rest frame) 
assuming a straight line for fq{X) (vertical axis) and a power-law 
(horizo ntal axis). The diagonal line is the line of equality, (see 
section 4.2). 


Another possible source of systematic error is in our se¬ 
lection of the wavelengths Xmax and Xmin, within which the 
suspected 4000A break and our figure of merit (Eqn. ^ is 
defined (see Fig. As discussed in section 3.1, the values 
Xmin — 3500A and A max = 5080A were chosen as a com¬ 
promise between: hrst, to maximise the number of sources in 
which the 4000A break remains within the observed wave¬ 
length after redshifting, and second, to make optimal use of 
the break region. What are the effects on the galaxy fraction 
if a different wavelength range were assumed? 

To explore this, we choose to vary Xmax alone. Due to 
the relatively small wavelength range at A < lOOOA avail¬ 
able in the galaxy spectrum (Fig. |^, we are not as flexible 
in varying Xmin- We therefore keep Xmin hxed at 3500A. 
Assuming the same bin widths (200A and 400A ) centered 
on Xmin and Xmax, and the PL dehnition for fq{X), we have 
computed galaxy fractions at 5000A with Xmax = 5500A , 
and Xmax = 6500A . Combined with our estimates using 
Xmax = 5080A , we hnd that changing Xmax makes neglible 
difference in the galaxy fraction. The fractions differ by no 
more than 1%. 

We now compare these systematic uncertainties with 
estimates of the random errors. Fig. ^ shows the distribution 
of random errors in the fractional galactic contributi on a t 
5000A as determined from our algorithm (see section |2.3| ). 
The range in random errors is signihcantly broad, with a 
majority of values exceeding our maximum Ict systematic 
uncertainty of 0.02 deduced from Fig. In the remaining 
sections, we therefore quote all uncertainties in the relative 
galaxy contribution as purely statistical, based on random 
errors alone. 


4.3 Galactic Contribution to the Observed 
Optical-Near IR Continuum. 

Using our algorithm, we have computed the fractional 
galaxy contribution in the individual bandpasses Bj (A ~ 
4400A) and A (A ~ 2.2/rm), expected in an observer’s frame 
for each source. We estimate these using Eqn. where huxes 
are approximated by computing the median spectral flux at 
the central wavelength of each bandpass. Since our observed 
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b0131m001 z = 0.8790 


b0338m214 z = 0.0480 




Rest Wavelength (A) 


Rest Wavelength (A) 


b0153m410 z = 0.2260 


b0406m31 1 z = 0.0565 




Rest Wavelength (A) 


Rest Wavelength (A) 


b0221p067 z = 0.5100 


b0454p066 z = 0.4050 




Rest Wavelength (A) 


Rest Wavelength (A) 


Figure 5. Observed rest frame spectra (dotted) for a number of sources in the Parkes sample. Also shown are reconstructed model spectra 
(solid) as defined by Eqn. ^ and power-law fits representing the underlying quasar continuum (dashed). Since this quasar spectrum is 
purely defined within the range 3500A(A^i„) < A < 5080A(Amaa:), the model spectra are only effective in reproducing those regions 
observed within this wavelength range. 


spectra do not extend to near-IR wavelengths, we estimate 
the observed spectral flux in K using our observed Bj — K 
colours and extrapolate from the spectral flux corresponding 
to Bj. Unfortunately, not all of the 53 sources used in our 
algorithm have a measured K magnitude and hence Bj — K 
colour. Of these sources, exactly 34 have known K magni¬ 
tudes. 

The distributions in galaxy fractions in Bj and K are 


shown in Figs. ^ and b respectively. Sources with mean 
galaxy fractions <5% are replaced by their 3a upper limits 
(dashed histograms). The distributions in Figs. and b 
appear very similar, except however for a greater number 
of sources with galaxy contributions > 70% in the 7F-band. 
These are all low redshift sources with strong 4000A breaks 
in their spectra. Their light is expected to be dominated by 
evolved stellar populations, and hence strongest in K. We 


© 1998 RAS, MNRAS 000, 000-000 






























8 F. J. Masci et al. 


Random Errors 



o(gal frac. at 5000 Ang.) 


Figure 7. Distribution of Itr random errors in the fractional 
galaxy cont ribution at SOOOA (rest frame) for Parkes sources, (see 
section |4.2|). 



Galaxy frac. in Bj 


must also note that the Bj magnitudes used to estimate the 
K sp ectral fluxes via Bj — K colours are only accurate to ~ 1 
mag (Drinkwater et al. 1997). This photometric uncertainty 
is likely to contribute significant scatter in our estimates of 
the K galaxy fraction in Fig. |^. 

We have divided the Tf-band identifications into two 
populations: those which show extended (resolved) struc¬ 
ture, and those which remain unresolved. The galaxy frac¬ 
tion in K for these two populations is shown as a func¬ 
tion of g in Fig. Open symbols represent resolved sources, 
and filled symbols, unresolved sources. As expected, those 
sources exhibiting resolved structure are also those which 
show large galaxy fractions and are at relatively low red- 
shifts. 

We now investigate whether emission from the host 
galaxies of Parkes quasars can significantly contribute to 
their observed Bj — K colours. We do this by computing 
the Bj — K colour of the hypothesised underlying “quasar”, 
{Bj — K)q, we would expect if contribution from the host 
galaxy was absent in each source. If the observed colours 
were entirely due to galactic emission, then we expect the 
distribution in {Bj — K)q to show a relatively small scatter, 
i.e. similar to that observed for optically-selected quasars 
where typically {Bj — K)q ~ 2.5. 

The colour of an underlying quasar, {Bj — K)q, can 
be written in terms of the observed colour {Bj — K)obs and 
the galaxy fractional contributions Fgai{Bj) and Fgai{K) as 
follows: 


(Bj - K)q 


{Bj - K)oba + 2.5 log 


^-Fgal{K) 
1 - Fgal{Bj) 


(13) 


{Bj — K)q is plotted against {Bj — K)obs in Fig. p^ . As can 
be seen, the scatter in galaxy subtracted colours, {Bj — K)q, 
remains and is extremely similar to that of the observed 
colour distribution. We quantify the galaxy contribution to 
the observed Bj — K colours (in magnitudes) in Fig. p^. 
From Figs. §a and b, we conclude that the observed spread 
in colours cannot be due to emission from the host galax¬ 
ies of Parkes quasars. An independent mechanism must be 
involved. 



Galaxy frac. In K 


Figure 8. (a) Distribution in fractional galaxy contributions in 
Bj and (b) K for 2 < 1. Dashed portions represent Sit upper 
limits, (see section |4.3|). 
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Figure 9. Fractional galaxy contribution in iC as a function of 
redshift for resolved sources (extended on K and Bj images; open 
symbols) and unresolved sources (closed symbols). Triangles rep¬ 
resent 3cr upper limits on the galaxy fraction. 
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Gal. contribution to observed Bj-K (mag.) 


Figure 10. (a) Observed Bj — K colour vs. colour of the underly¬ 
ing quasar (i.e. resulting galaxy subtracted colour: (Bj—K)q). Er¬ 
ror bars correspond to Itr significance, (b) Distribution in galaxy 
contribution to observed Bj — K colours in magnitudes. 


5 A TEST FOR THE UNIFIED MODEL 


Studies of the host galaxy properties of BL-Lacs and radio¬ 
quasars can be used as a test of the unified scheme for radio- 
loud AGN. Motivated by the canonical axisymmetric model 
for AGN, the basis of this scheme is that the appearence 
of an extragalactic radio source is primarily determined by 
viewing geometry. Extended FRI and FRII-type radio galax¬ 
ies are believed to represent the parent (misaligned) popu¬ 
lations of t he more compact BL-Lacs and radio quasars re - 
spectively (Urry & Padovani 1995 and references therein). 


If classification is purely based on orientation, then intrinsic 
properties such as host galaxy luminosity should be approx¬ 
imately uniform throughout. In this section, we shall test 
this hypothesis. 

Studies have shown that for quasars and FRII radio 
galaxies at redshifts a<0.3, the situation is not entirely clear. 
From comparisons of their mean host galaxy luminosities, 
some studies have shown that FRII h osts are fainter by 
~ 0.5-1 mag. (smith & Heckman 1989), while others have 


concluded that they are comparable (eg. Taylor et al. 1996). 
The main difficulty in these studies was finding sufficiently 
large samples of radio galaxies and quasars matched both 
in radio power and redshift. There is strong observational 


evidence however that the low redshift BL-Lacs reside in 
giant ellipticals with mean optical luminosities and d e Vau- 
coleurs law profiles similar to those in FRIs (Ulrich 


1988; Stichel et al. 1993). Very little is known about the 


host galaxies of compact radio sources at higher redshifts. 
Using our algorithm however, we can get estimates of host 
galaxy K magnitudes for sources up to « ~ 1. 

We can predict the host galaxy K magnitudes of Parkes 
sources directly from our estimates of the fractio nal g alaxy 
contribution in the A-band, Fgai{K) (see section 4.3). The 


host galaxy magnitude in an observer’s frame can be written: 


^source 


- 2.5\og[FgaiiK)], 


(14) 


where Kaourcs is the o bser ved K magnitude of the source. 
As discussed in section 4.3, Fgai(K) is determined from the 
median spectral flux at A ~ 2.2/rm (Eqn. which we es¬ 
timate using the observed spectral flux at A ~ 4400A and 
Bj — K colour. 

Estimates of Kgai as a function of redshift are shown in 
Fig. 0- For comparison, we also show the range observed for 
radio galaxies (shaded region) as determined from a num¬ 
ber of independent studies (McGarthy 1993). Within our 
quoted uncertainties, there appears to be no significant dif¬ 
ference in the mean host galaxy magnitude of “compact” 
Parkes sources and extended radio galaxies at the redshifts 
indicated. The compact sources however appear to show a 
larger scatter in Kgai at some redshift. Since the Kgai val¬ 
ues were determined from non-contemporaneous measures 
of observed Bj spectral fluxes and Bj — K colours, this may 
be attributed to variability in the underlying AGN. We are 
unable at present to quantify this uncertainty. From Fig. 
we conclude that the host galaxy luminosities of these two 
classes of radio source is consistent with that required by 
the unified model. 


6 DISCUSSION 


Our results of section 4.3 clearly show that galactic emis¬ 
sion is unlikely to fully explain the dispersion in Bj — K 
colours observed. This conclusion is only valid however for 
sources at redshifts 2<0.9. At higher redshifts, the 4000A 
break feature on which our algorithm is based is redshifted 
out of our observational wavelength range. Taking into ac¬ 
count our completeness in spectral data (section only 
53 of the 323 sources in the Drinkwater et al. (1997) sam¬ 
ple have been analysed using our algorithm. Since we are 
limited to redshifts 2:;<0.9, it is possible that we are biased 
towards detecting relatively large galaxy contributions. Our 
results may thus not be representative for the whole sample 
of Parkes quasars. 

We have two strong pieces of observational evidence 
that supports a minimal galaxy contribution from the high 
redshift Parkes sources: First, all sources with 2>0.5 appear 
very compact in K (eg. Fig. and exhibit broad-line equiv¬ 
alent widths typical of those observed in optically-selected 
quasars. Second, significantly high levels of linear polari¬ 
sation (^5%) have been observed in the near-IR in several 
sources at 2:>1 (see Masci 1997). This strongly indicates that 
the emission is dominated by a non-thermal mechanism. It 
is important to note however that the number of sources 
in which a polarisation has been searched for is too low to 


© 1998 RAS, MNRAS 000, 000-000 














































10 F. J. Masci et al. 
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Figure 11. Host galaxy K magnitude as a function of redshift 
for Parkes sources (symbols) and the range observed for radio 
galaxies (shaded region; McCarthy 1993). 


draw any reasonable conclusion. Further polarimetric stud¬ 
ies of preferably the reddest quasars are necessary to asses 
the importance of a host galaxy component. 

The next significant step in improving the algorithm 
presented here for determining the galaxy contribution 
would be to completely discard our assumption of a power- 
law for the shape of the underlying quasar spectrum. Al¬ 
though a simple power-law is a good overall representation 
for the continuum around A ~ 4000A, generic quasar spectra 
also include a complex blend of emission line features super¬ 
imposed on the continuum about this region. As seen in most 
compilations of composite QSO spectra (eg. Francis et al. 
1991), wavelengths shortwards of 3700A are contaminated 
by a Balmer coninuum, where emission from the convergence 
of high order Balmer lines can introduce a steep rise in this 
part of the spectrum. Furthermore there is also a complex 
blend of Fell emission setting in at these wavelengths. These 
considerations therefore invalidate our assumption of a pure 
power-law for the “smooth” quasar spectrum. Nonetheless, 
it is likely that the omission of these additional intrinsic 
features at A<4000A will have led us to overestimate the 
galaxy contribution on average. 


7 CONCLUSIONS 

This paper has explored whether emission from the host 
galaxies of Parkes quasars can significantly contribute to the 
relatively large spread in Bj — K colours observed. If the 
hosts are classical giant ellipticals and their flux strongly 
contributes, then this would be expected since elliptical 
colours are known to be quite red in B — A to 2 ~ 2. 

We have devised an algorithm that measures the rela¬ 
tive galaxy contribution in each source in an unbiased way 
using the characteristic 4000A break feature of elliptical 
galaxy SEDs. The basis of the algorithm involves subtract¬ 
ing a generic elliptical SED from each source spectrum un¬ 
til the 4000A break feature disapears and what is left is 
a “smooth” spectrum containing no breaks. This “smooth” 
spectrum we refer to as the underlying quasar spectrum. 
The only requirement by our algorithm is that this remain¬ 
ing spectrum be smooth. The galactic contribution, relative 


to the total light at any wavelength is estimated from the 
amount of galaxy subtracted. 

The main conclusions are: 

1. For z^O.9, (for which the 4000A feature remains ob¬ 
servable in our spectra), we Hnd broad and almost bimodal 
distributions in the relative galaxy fraction in Bj and K. 
Most sources (^70%) have galaxy fractions < 0.3 at the 3a 
level in both Bj and K. The remainder have large galaxy 
contributions and are predominately low redshift galaxies 
with strong 4000A breaks. All of these latter sources are 
spatially extended and resolved on Bj and A-band images. 
In particular, there is a clear distinction in the strength of 
the 4000A break for resolved and unresolved sources. 

2. Using these estimates, we End that the mean A- 
band magnitude of the host galaxies of flat spectrum radio 
quasars is consistent with that of extended radio galaxies at 
z^O.9. This is consistent with the unified model for radio- 
loud AGN. 

3. By subtracting the galaxy contribution in each band¬ 
pass from the observed Bj — A colours of Parkes sources, 
we End that at predominately the 2a confidence level, the 
relatively large spread in colours still remains. We conclude 
that in a majority of cases, the relatively red colours must 
be due to a mechanism other than that contributed by a 
“red” stellar component. 
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